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ABSTRACT 


A  new  method  has  been  developed  to  determine  molecular  weight  (specifically,  ratio 
2 

g  =  soj  /D,  where  ou  =  angular  velocity,  s  =  sedimentation  coefficient,  D  =  diffusion  coeffi¬ 
cient)  from  an  experimental  schlieren  plot,  dc/dr  vs  r,  taken  at  equilibrium.  Computation  of 
the  quantity  g  can  be  accomplished  by  applying  the  newly  derived  formula: 

(  -jr  -^r  )  =  irv( gk)  +  -g-  gr2, 

where 

k  =  (c )  exp  ( — ^-gr2) 


is  a  constant,  and  c  is  concentration  at  distance  r  from  center  of  rotation. 


For  a  case  of  apolydisperse  system,  the  following  method  of  computation  has  been  developed. 
The  heterogeneous  system  was  replaced  by  an  equivalent  N-component  paucidisperse  solute. 

A  set  of  N  values  was  arbitrarily  chosen,  where  each  g^  is  related  to  component,  i, 
present  in  the  system.  The  following  expressions  were  defined: 


y.  =  2  /  r .  (  r  2  - 

I  J  b 


is  a  quantity  measured  at  point  j  on  the  experimental  plot;  b  and  m  denote  bottom  and  menis¬ 
cus,  respectively;  and  the  function 


a  --[a1'1 


Since  it  is  shown  that 


where  cq^  is  the  initial  concentration  of  component  i,  evaluation  of  a„  and  y^  experimentally 

leads  directly  to  the  determination  of  c  Therefore,  by  applying  this  method,  one  can 

determine  the  entire  distribution  of  molecular  weight  from  equilibrium  sedimentation,  as 
well  as  the  average  values. 
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SECTION  I 
INTRODUCTION 


1.  EQUILIBRIUM  SEDIMENTATION 

Equilibrium  sedimentation  was  first  applied  by  Svedberg  in  his  early  works  (Reference  1). 
Since,  at  equilibrium,  concentration  of  polymers  in  solution  does  not  depend  upon  time,  the 
equation  of  ultracentrifuge  reduces  to  an  ordinary  differential  equation: 

dc 

- —  =  -  div  J  =  0  ,  ( I  ) 

dt  — 

where 

c  =  polymer  concentration 
£  =  flux  of  solute. 

Hence 

|jJ  =  I  -  Dcr  +  stn2cr  I  =  const.  (2) 

D  and  s  are  coefficients  of  diffusion  and  sedimentation,  respectively. 

cj  =  angular  velocity  of  ultracentrifuge; 

r  =  distance  from  center  of  rotation 

c  =  dc/dr. 
r 

Since  at  equilibrium  no  flux  occurs,  jJ  =  0,  the  following  formula  holds: 

.  dc  2  ,,, 

D  -t—  =  sco  cr  (3) 

a  r 

2 

By  introducing  a  new  constant,  g  =  sen  /D,  the  Differential  Equation  (3)  can  be  written  in  a 
short  form: 


dc 

dr 


=  gre, 


(4) 


or; 

d(^i/c)=  —  gcd(r2)  (5) 

2 

Formula  (5)  leads  to  a  well  known  method  of  computation  (Reference  1)  which  determines  g 
from  any  two  points  i  and  j  of  an  experimental  plot  c  vs  r: 

c./c.)  =  -^  g  (  r.2-r  2  ).  (6) 
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By  substituting  g  into  the  Svedberg  Equation  (Reference  2),  molecular  weight  can  be  easily 
evaluated: 


M 


RTg 

(I  -Mp  )tu2 


(7) 


where 

M  =  molecular  weight 
R  =  gas  constant 
T  =  absolute  temperature  in  °K 
V  =  partial  specific  volume  of  solute 
p  =  density  of  solvent. 

2.  LIMITATIONS  OF  EQUILIBRIUM  SEDIMENTATION 

The  following  are  limiting  factors  in  the  application  of  equilibrium  sedimentation. 

a.  The  time  of  the  experiment  becomes  exceedingly  long  if  M  >  250,000.  If  M  >  1,000,000 
this  method  is  not  feasible.  Therefore,  equilibrium  sedimentation  can  only  be  advantageously 
applied  if  molecular  weight  is  below  250,000. 

b.  Schlieren  optics  is  the  simplest  and  most  generally  used  experimental  technique.  How¬ 
ever,  this  method  provides  a  plot  ofdc/drvs  r.  Since  the  conventional  methods  of  computation 
require  a  plot  of  c  vs  r  (see  Equation  (6)),  the  experimental  schlieren  curve  has  to  be  trans¬ 
formed  into  an  integral  one.  Such  a  transformation  is  very  laborious  and  not  sufficiently 
precise. 

c.  A  plot  of  c  vs  r  can  be  obtained  by  applying  interference  optics.  Unfortunately,  this 
method  cannot  be  used  where  high  speed,  elevated  temperature,  or  corrosive  solvents  are 
involved. 

d.  Absorption  optics  also  provide  a  plot  c  vsr.  However,  this  technique  requires  a  tedious 
procedure  of  calibration  and  evaluation,  and  very  good  conditions  for  precise  photographic 
processing. 

e.  In  recent  years  a  new  technique  has  been  under  development  (Reference  3)  which  also 
leads  to  an  experimental  plot  of  c  vs  r.  This  method  is  based  upon  direct  scanning  of  an 
ultracentrifugal  cell  by  photosensitive  elements.  Unfortunately,  this  optical  system  is  not  yet 
complete  and,  in  addition,  necessary  equipment  is  very  expensive. 

3.  PRIOR  ATTEMPTS  TO  IMPROVE  EQUILIBRIUM  SEDIMENTATION 

Some  attempts  have  been  made  to  improve  the  method  of  equilibrium  sedimentation  and  to 
relax  its  limitations. 

Archibald  (References  4  and  5)  applied  Equation  (4)  to  evaluate  g  at  the  meniscus  at  any 
time;  e.g.,  before  equilibrium.  Since  no  transport  can  exist  at  the  meniscus  and  bottom  of  the 
ultracentrifugal  cell,  an  equilibrium  equation  must  always  hold  at  this  particular  point.  Some 
of  his  evaluations  are  of  general  use,  and  are  mentioned  subsequently.  Experimental  difficul¬ 
ties,  however,  arose  which  limited  the  application  of  Archibald’s  method. 
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Another  method,  “sedimentation  at  constant  density  gradient,”  was  successfully  developed 
by  Meselson,  et  al.  (Reference  6).  This  technique,  which  significantly  shortens  the  time  of  an 
experiment,  requires  the  presence  of  a  special  binary  or  multicomponent  solvent.  Such  a 
system  is  supposed  to  give  a  constant  density  gradient  and  to  cause  sedimentation  and  flo¬ 
tation  of  polymer  toward  the  middle  of  the  cell.  It  is  evident  that  each  family  of  polymers 
can  be  investigated  in  this  manner  only  if  an  appropriate  solvent  system  has  already  been 
developed.  Unfortunately,  at  present  there  are  only  two  successful  systems  —  one  for  deoxy¬ 
ribonucleic  acid  and  one  for  polystyrene. 

Since  molecular  weights  of  polymers  under  current  investigation  in  the  Air  Force  Materials 
Laboratory  do  not  exceed  250,000,  conventional  equilibrium  sedimentation  can  be  successfully 
applied.  Therefore,  the  objective  of  this  work  is  to  improve  the  methods  of  computation  and 
determine  g  directly  from  an  experimental  plot  of  cr  vs  r. 


SECTION  II 

MATHEMATICAL  ANALYSIS  OF  EQUILIBRIUM  SEDIMENTATION 


Integration  of  Differential  Equation  (4)  leads  to  the  following  solution: 

2 


c  =  k 


2  gr 


(8) 


where  k  is  a  constant 
Hence: 


i  2 


i  2 


”Tgf2  “Tgrm  “Tgrb 

k  =  ce2  =  cme  m=  cbe 


(9) 


m  denotes  meniscus  and  b  bottom  of  the  cell. 

Following  Equation  (9)  the  concentration  at  any  point  can  be  given  as  a  function  of  any  chosen 


reference  point  (c^,  r^) : 


if2  2  \ 

c  =  c  e  2  H 


(10) 


The  first  derivative  of  c  is  given  by  Equation  (11): 


cr  =  grcq  Q 


i  i  2  2\ 

T  Q(  r  ~fq  ) 


(ID 
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Further  differentiation  of  Equation  (11)  leads  to  the  following  relationship: 

I  ,2  2.  I  .  2  2  , 

-g  g(  r  -  rq  )  2  2  T  O'  r  ~  rq  ' 

crr=  gcqff  +  g  r  cq«? 

By  combining  Equations  (10),  (11),  and  (12),  one  will  obtain: 


or: 


crr=  (  T  +  9r)cr 

d  jQsn/  Cr  I 


d  r 


=  g' 


Equations  (8)  and  (14)  lead  to  the  following  expression: 


(12) 


(13) 


(14) 


-£w  cr  =  ^/rv  (kgr)  +  gr2  (15) 

and  finally: 

^rv  (Cr/r  )  =  in/  (gk)  +  y  gr2.  (16) 

A  2 

Equation  (16)  can  be  applied  for  the  computation  of  g.  By  plotting  J Kru  (c^/r)  vs  r  ,  one  can 

obtain  the  constant  values  g/2  and  Hrc(gk) .  With  the  help  of  g  and  k  the  concentration  c  can  be 
evaluated  at  any  r;  cm  and  c^  (concentrations  at  meniscus  and  bottom  respectively)  can  be 

determined  in  the  same  manner. 


There  is  another  way  to  evaluate  g  by  approximation  in  the  vicinity  of  a  chosen  reference 
point  (cf 

tion  (10)  reduces  to  the  following  series: 


2  2 

_  ,  r  ).  Since  under  these  circumstances  1/2  g(r  -  r  )  is  a  small  number,  Equa- 

q  q  q 


(17) 


2  2  2 
A  plot  of  (c  /r)  vs  (1/2)  (r  -  r  )  provides  values  of  the  constants  g,  c  ,  and  g  c  ,  and  makes 

it  possible  to  determine  g  and  c  .  Since  no  flux  exists  at  the  boundaries  of  an  ultracentrifugal 

q 


cell,  any  equation  describing  equilibrium  always  holds  at  the  meniscus.  Therefore,  the  method 
of  computation  expressed  by  Equation  (16)  can  be  applied  in  the  vicinity  of  the  meniscus,  and 
such  an  application  can  be  made  at  any  time;  e.g.,  even  before  equilibrium  is  established. 
Under  these  conditions  g  (and  c  )  have  to  be  determined  as  limiting  values  for  r  — 
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SECTION  m 

POLYDISPERSE  IDEAL  SYSTEM 


Let  us  consider  an  ideal  dilute  polydisperse  system  composed  of  N-noninteracting  com¬ 
ponents.  At  equilibrium  each  component  attains  its  own  equilibrium  state.  Therefore,  there 
exists  a  set  of  the  following  N  equations: 


(i)  (i)  (i)  2  (') 

D  cr  =  s  o>  rc 


where  (i)  denotes  component  i. 

Or  in  terms  of  g^  =  ai  : 


(i)  (i)  (i) 

v  =  9  rc 


Solutions  of  the  differential  equations  given  by  Formula  (19)  can  be  expressed  either  by 
constants  or  by  any  reference  concentration  c^. 

,*%  i*\  i  ( i )  2  /.%  1  (i),  2  Zx 

k(l  r  =  c<q'eT5  <r  "'I1-  U  ,">)  =  120) 

Another  relationship  can  be  derived  from  the  law  of  conservation  of  matter.  Since  the  amount 
of  solute  enclosed  in  a  cell  remains  constant,  the  following  equation  holds  at  any  time: 


C  b 

=  J  'I'cdr, 


where  A  is  the  amount  of  solute  in  the  cell,  and  \f/  (r)  is  the  area  of  the  cylindrical  surface 
within  the  cell  at  the  distance  r  from  the  center  of  rotation.  In  the  case  of  a  sectorial  cell 
with  angle  <f>  and  depth  h: 

f  =  <f>  rh  (22) 

and  Equation  (21)  goes  into  the  following  formula: 


c  b 

/  h4> 


crdr 


where  A,  h  and  (f>  are  constant. 

If  c0  denotes  the  initial  homogeneous  concentration,  one  can  derive  from  Equation  (23)  the 
following  relationship: 


c  (  r3  -  r  2 )  =  J 

o  b  m  J . 
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Equations  (21  through  24)  were  first  derived  by  Archibald  (References  4  and  5). 
By  substituting  Equation  (8)  into  (24)  one  can  obtain: 


Hence: 


c 

o 


I  2  2 

c  (  r.  -  r 
o  d  m 


=  -f  Cch 
9  b 


c  ) 
m 


(25) 


(26) 


Equation  (26)  was  derived  in  a  different  way  by  Van  Holde  and  Baldwin  (Reference  9). 

The  value  (c^  -  c  )  which  expresses  the  increase  of  concentration  between  the  meniscus 

and  bottom  of  the  cell  can  be  easily  determined  from  experiment  by  integrating  c  between 
all  boundaries:  1 


c 


b 


c  = 
m 


cfdr  . 


m 


(27) 


By  substituting  the  value  of  (c^  -  cj  into  Equation  (26)  one  can  easily  determine  g. 

In  the  following,  Equation  (26)  will  be  applied  to  evaluate  the  distribution  of  a  polydisperse 
system.  Since  Equation  (26)  holds  for  any  component  of  a  heterogeneous  system,  one  can  write: 


J  <i)  (i>l 

(i)  2icb  -  cm  j 

°  ’  0(i>(rh2-r^) 

"  b  m 


(28) 


Equations  (20)  and  (11),  if  applied  at  the  meniscus  for  component  i,  go  into  the  following 
expressions: 


and: 


(i) 


g(i)r 


G 


(i).  2 
g  (r  - 


2 

rm 


(29) 


(30) 


First,  the  quantity  will  be  eliminated.  By  combining  Equations  (28)  and  (29)  the  following 
expression  will  be  obtained: 


(i) 


1  '  rb  rm  ' 


.  ti).  2 
T9  (  rK  ~  r 


m 


(i) 
1  m 


(31) 
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By  eliminating  the  quantity  c  ^  from  Equations  (30)  and  (31),  one  may  express  cr^  as  a 

function  of  g^  and  c  ^ : 

&  o 


c 


(i) 

r 


r  (  r 


2 

b 


2 


,(i)2*T3 


( i )  2 
r 


-jrg^r.2  4q(1V 


e 


e 


2  9' 


m 


(32) 


Formula  (32)  can  be  applied  to  any  one  of  the  N  components  present  in  the  system  under 
investigation.  Therefore,  this  formula  represents  a  set  of  N  equations.  By  carrying  out  the 
summation  of  all  N  equations  of  Formula  (32)  one  will  obtain  the  following  relationship: 


1  ,  2  2»  v 

T r  rb  m 1  £ 

i  =  I 


2  i  (i)  2 

g(i,%^9  r 


I  (i)  2 


_  T  9  rb  T  9 
Q  Q 


i  (i)  2 


m 


(i) 


(33) 


where  cr 


2,  cr  can  be  directly  measured  experimentally, 
i  =  I 


The  objective  of  this  section  is  to  evaluate  the  molecular  weight  distribution.  Since  the 
ratio,  P: 


P 


RT 

(  I  -  Wp)  o)  2 


(34) 


which  appears  in  Equation  (7)  is  an  experimental  constant,  identical  for  all  components 
present  in  the  system: 

(  i) 

M  =  Pg  ,  (35) 


the  distribution  of  molecular  weight,  M^,  can  be  replaced  by  a  distribution  of  quantities  g^. 
It  should  be  mentioned  that  evaluation  of  a  distribution  means  an  appropriate  adjustment  of  the 

pairs  (g^,  oo«).  Such  a  definition  gives  one  the  freedom  to  arbitrarily  choose  sets  of 

quantities  mentioned  above.  This  freedom  is  utilized  in  the  following  way: 

First,  a  set  of  N  g^  values  will  be  chosen  arbitrarily.  However,  one  has  to  bear  in  mind, 
that  the  experimental  average  g  values  cannot  be  outside  of  this  set.  In  addition  such  a  set 
cannot  be  very  narrow,  and  it  must  be  extended  over  a  wide  range  of  possible  g  values. 
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Let  us  now  denote  by  j  a  specifically  chosen  point  (r.,  c  )  on  the  experimental  curve.  For 

j  (A 

any  chosen  pair  (r.,  cr  )  and  for  any  arbitrarily  chosen  value  of  g  ,  the  following  quantities 
can  be  determined: 

^  .  i  /  ^  ^  \  Z-'  \ 


y.  =  2c  .  /  r .  (  r  -  r  ) 

’ j  ri  j  b  m  • 


a. .  = 
Ji 


.2  I  _  (  i  )  2 

9(,le?9  ri 

I  l  i  )  2  |  (i)  2  * 

T  9  fb  T9  r  m 


By  introducing  the  definitions  given  in  Equations  (36)  and  (37),  the  expression  cited  in  Equa¬ 
tion  (33)  can  be  transformed  into  summation: 


where  y^  and  a„  are  known  quantities. 

Let  us  choose  L  pairs  (r.,  c  )  where  L>  N.  Since  Equation  (38)  holds  for  any  one  of  such 

3  Pj 

defined  L  pairs,  there  exists  a  system  of  L  simultaneous  linear  equations,  with  N  unknown 
quantities  c 


y'  =  *«  c° 


»L  *  I.  °Li  cc 

I  =  I 


The  following  definitions  of  two  vectors  and  one  matrix: 


yL>  ; 


co  3  • 


(4  2) 
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make  it  possible  to  express  the  System  indicated  in  Equation  (39)  in  a  compact  form: 

_Y  =  AC  .  (43) 

If  L  =  N  one  can  immediately  write  the  solution 

C_  =  B_Y  .  (44) 

where  B  =  A  *  is  the  inverse  matrix  of  A. 

If  L  >  N  a  method  of  least  squares  can  be  applied  to  evaluate  vector  C. 

Since  the  elements  of  C  are  correlated  to  the  appropriate  set  of  values,  the  evaluation 
of  vector  _C  accomplishes  the  determination  of  the  molecular  weight  distribution. 

The  theory  derived  in  this  work  is  an  approximate  one.  It  disregards  any  influence  of  virial 
coefficient,  interactions  within  a  multicomponent  system,  compressibility  in  the  ultracentrifuge, 
and  other  nonspecified  factors. 
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is  a  quantity  measured  at  point  j  on  the  experimental  plot;  b  and  m  denote  bottom  and 
meniscus,  respectively;  and  the  function 
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Since  it  is  shown  that 


y.  = 
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where  cQ'^  is  the  initial  concentration  of  component  i,  evaluation  of  a.,  and  y.  experimen¬ 
tally  leads  directly  to  the  determination  of  c  Therefore,  by  applying  this  method,  one 
can  determine  the  entire  distribution  of  molecular  weight  from  equilibrium  sedimentation, 
as  well  as  the  average  values. 


